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Abstract. Intense light with frequencies above typical atomic or molecular ionization 
potentials as provided by free-electron lasers couples many photons into extended 
targets such as clusters and biomolecules. This implies, in contrast to traditional 
multi-photon ionization, multiple single-photon absorption. Thereby, many electrons 
are removed from their bound states and either released or trapped if the target charge 
has become sufficiently large. We develop a simple model for this photo activation to 
study electron migration and interaction. It satisfies scaling relations which help to 
relate quite different scenarios. To understand this type of mult i- electron dynamics on 
very short time scales is vital for assessing the radiation damage inflicted by that type 
of radiation and to pave the way for coherent diffraction imaging of single molecules. 
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1. Introduction 

The difference between multi-photon ionization at near-infrared frequencies (~1.5eV) 
and at extreme-ultraviolet (XUV) frequencies or higher (>50eV) lies in the fact that 
multi-photon sing/e-electron ionization is drastically suppressed in the latter case due to 
the small dipole matrix elements for continuum-continuum transitions which typically 
would be required for an electron to absorb more than one photon. As a consequence, 
multi-photon ionization by XUV or X-ray radiation, as available from novel free-electron 
laser sources [1-3], means smg/e-photon ionization of many atoms in an extended target 
such as a cluster or bio-molecule [4-6]. We call this process photo activation. As other 
absorption mechanisms, such as inverse bremsstrahlung, are less important for high laser 
frequencies, the system is "driven" by this activation process only. In the following we 
will develop a fairly general and simple model, which we call Coulomb complex (CC). It 
focuses on the multi-electron dynamics treating the ions created as a homogeneous and 
static background charge. The model shares similarities with the shell model for nuclei 
or metal clusters [7] with the difference that the electrons are through photo activation 
far above the ground-state which allows us to follow their dynamics classically. 

We will first discuss the initial multi-electron ground state in Sect. 2. The ground 
state is the prerequisite for photo activation which we specify in Sect. 3 including the set- 
up of a universal time-dependent photo activation rate. In Sect. 4 we discuss a powerful 
scaling property of the photo-activated CC which allows one to relate quite different 
situations in terms of photo activation time (length of the XUV pulse) T, excess energy 
e* of the electrons after photo absorption and size of the target designated by a radius 
R, to each other. In Sect. 5 we elaborate on a generic example and discuss relations of 
the present model to recently measured electron spectra of Xenon clusters [5] exposed 
to 90 eV photon energy pulses at FLASH [1]. In Appendix A we give an analytical 
expression for the electron spectrum for the case of sequential ionization. 

2. The Coulomb complex before photo activation 

We assume a system with N electrons which will be eventually photo-activated. 
They are bound by positive ions which are assumed to be fixed in space during the 
relevant electron dynamics. Moreover, we approximate the ions within the CC as a 
homogeneously charged sphere, with total charge Q and radius R, producing a radial 
potential for the electrons of the form 



The potential is within the CC (r < R) harmonic with the frequency = (Q/R 3 ) 1 / 2 and 
for r > R of pure Coulomb nature. The full Hamiltonian of the interacting electronic 
system is given by 




for r < R 
for r > R. 



(1) 




(2) 
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with rjk = \rj — Tk\ and the prime excluding the term j = k. The total ground state 
energy E gs may be approximated by considering a homogeneous electron distribution. 
Then we can replace the sum over all electrons by an average potential energy v c and 
an average electron-electron interaction energy v ee , i.e., 

n 1^1 N 

E gs = E Hrj) + 2 E'— = Nv c + -v ee . (3) 
j=i j,k=i i* 

This assumption should be valid for sufficiently many electrons N, but applies also for 
quite small iV as we will see below. Within this approximation the average energy from 
the background potential reads 

«c=^//W) = -^. (4) 

The interaction energy v ee cannot be calculated this way, since we have to exclude 
explicitly the self-interaction term. Instead we estimate the average energy of an electron 
e = v c + v ee in the ground state. This energy can be approximated under the assumption 
that the electron density is identical to the ion density apart from a "hole" with volume 
(4:Ti/3)R 3 /N. The corresponding potential energy of such a hole is most easily calculated 
for a sphere (with a radius R/N 1 / 3 ) and gives 

3iVl/3 ^ 

This single-particle binding energy depends exclusively on the density of the CC. Finally 
we can calculate the average electron-electron repulsion energy 

3 7VV3 6Q 

v - = e - v ° = -2nr + 5R> (6) 

which leads to the explicit approximation 

N 3N 2 3 iV 4//3 

-C/ es = NV C H V ee = , 7 

gs 2 ee 5 R 4 R , \ ) 

for the total energy of Eq. (3) of a neutral CC with N = Q. Note, that due to the non- 
linear dependence of the average electron-electron interaction upon the electron density, 
.Egs does not scale quadratically with N. Table 1 presents the analytically estimated 
energies e and E gs from Eqs. (5) and (7), respectively, in comparison to the numerically 
obtained values for 9 different values of N which only for the smallest value N = 7 differ 
by more than 3 %. 

The numerical energies were obtained by propagating N electrons with the 
Hamiltonian (2) while reducing the particle velocities by a factor of 0.75 every 2.4 fs 
for N < 123 and every 0.24 fs for iV > 123 to arrive at a minimal energy configuration 
of the system with corresponding optimized positions {f j} for the electrons. The binding 
energies of the N electrons are given by 

Si = V(u) + E'^ . (8) 
fc=i ' lk 

Figure 1 shows these binding energies for various CCs with the same radius R but 
different electron numbers N. Interestingly, not only the average values e agree quite 
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Figure 1. Individual electron energies £j according to Eq. (8) as a function of the 
electron's distance r, = |r<| from the origin of the Coulomb complex for i? = 26 ao 
(Bohr radius) and various electron numbers N. The corresponding approximated 
average energies e from Eq. (5) are shown with lines. 

well with the analytical estimates (shown by lines in Fig. 1), there is also relatively small 
variation in the individual numerical electron energies (shown by symbols in Fig. 1). This 
applies down to such small electron numbers as iV = 7. In the figure one also clearly 
recognizes the formation of electron shells as known from so called Coulomb crystals [8] . 

3. Photo activation of the Coulomb complex 

In the previous section we have characterized the electrons of the CC in energy and 
prepared them in space to absorb photons. Due to the small dipole matrix elements in 
the continuum it is much more likely that N p h photons of XUV energy from an intense 
pulse are absorbed by N ph different bound electrons compared to the situation of multi- 
photon absorption by a few electrons. Hence, we have the N p ^ photons absorbed each by 
one bound electron. The single-photon absorption rate is proportional to the intensity 



N 


7 19 27 57 123 253 461 949 1935 


—e [eV] (num.) 

(anl.) 

—E gs /N [eV] (num.) 

(anl.) 


2.9 4.1 4.6 5.9 7.7 9.8 12.0 15.3 19.4 

3.0 4.2 4.7 6.0 7.8 9.9 12.1 15.3 19.6 

6.1 14.4 19.7 39.4 81.9 163.5 296.7 605.1 1226.7 
5.9 14.0 19.3 38.8 81.1 163.8 295.5 603.6 1224.9 



Table 1. Comparison of numerical (num.) and analytical (anl.), cf. Eqs. (5) and (7), 
binding (e) and ground state (E gs ) energies of a Coulomb complex with N electrons 
relaxed in a background potential of Eq. (1) with Q = N and radius R = 26 ao- 



: 
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I(t) of the pulse, 

^ = N a JM (9) 
dt u> 

with n(t) the number of electrons that have absorbed a photon up to time t, while 
iV at is the number of atoms and the photo-absorption cross section at frequency oo. 
Assuming for convenience (in fact one can use any pulse shape) a Gaussian pulse with 
full width at half maximum T we get 



n{t) = N ^ I ^f dt' exp(— 4 In 2(t'/T) 2 ) 



= ^[l + erf(2v^2|)], (10) 

with iVp h = NbxVujIqT (7r/4m2) 1//2 joj the total number of photons absorbed from the 
pulse. The expression for N p ^ requires that each atom is singly ionized by the pulse 
with the same cross section er^. This is of course an idealization, since realistically 
the cross section changes even if each atom is only singly ionized through ionization 
into a continuum which already contains electrons. Moreover, not the highest occupied 
molecular orbital may get predominantly ionized leading to more than one electron per 
photon through auto-ionization, and finally, simply more than one photon could be 
absorbed by each atom leading to multiple ionization. Nevertheless, Eq. (10) provides a 
reasonably general yet realistic form for the number of absorbed photons. 

Within our model we account for the photo absorption by the activation of electrons, 
which means that only after a certain time do the electrons take part in the dynamics 
with an initial momentum specified below. This implies that each electron j is held 
fixed at its original position until its individual time tj of activation. Since the photo 
absorption is a statistical process we treat the activation process statistically. Therefore 
we first calculate N p h random numbers {xj} between and 1. The activation time tj of 
electron j is then given by solving the equation n(t) = Xj iV p h which is implicit in t. In 
order to keep the computational expense as low as possible the number of electrons iV 
of the CC is set equal to the number of activated electrons: N = 7Vp h . For simplicity, 
we consider a CC which is neutral in the beginning, so that finally Q = N = N ph . 

We determine the initial momenta pj of the photo-activated electrons under the 
condition that activation of a single electron should lead asymptotically (if this electron 
is removed from the CC) to the atomic excess energy e*. This implies for the initial 
kinetic energy for each activated electron in the complex 

P 2 N 1 

| = e* - Vfr) - £ '- . (11) 

z fc=l ' jk 

For large N we can estimate the initial average kinetic energy per electron £km by 
averaging Eq. (11) over all electrons and obtain with Eq. (5) 

1 N 3 iV 1 / 3 

e ki n = ^ E Pj = £ * ~ fa + ^cc) = e* + - — . (12) 
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Since the numerical ground state energies for the individual electrons are quite 
homogeneously distributed this also carries over to the kinetic energies p^/2 which 
are close to the average values £kin- Note that u appears here only indirectly through 
determining e* = u — E- ip with the true atomic ionization potential E- ip . 

Formally we can express the activation by means of an activation function! 

A T (X,Y,t) = e(r-t)X + e(t-r)Y (13) 

in order to rewrite the original Hamiltonian (2) as 

H A =Y.\ A,0 % Pi ' t) \ v(A tl {f i ,r i ,t))j 

N i 




E'i • ( 14 ) 

2 j>k=1 1 A t . (Tj,Tj,t)-A tk (T k ,T k ,t) 

For this Hamiltonian the position of electron j in phase space is kept constant at (fj, pj) 
until its activation at time t — tj. One should note that, however, the single-electron 
energy Ej may change already for times t < tj due to the interaction with previously 
activated electrons. This is quite important since in highly charged systems the effects 
of the surrounding charges (binding, screening, etc.) easily exceed atomic properties by 
far and may change on a femtosecond time scale. The calculation of all interactions 
(including those with the not yet activated electrons) is therefore crucial. 

Note that, although the Hamiltonian Ha is time-dependent through the activation 
functions A, the corresponding total energy E(t) is conserved if one starts the 
propagation at the initial point in phase space, i.e., Tj(tj) = fj and Pj(tj) = pj for 
j = 1, . . . , N. It is lim^o [E(tj+6t) - E(tj-6t)] = for any j. 

4. Scaling in the dynamics of Coulomb complexes 

The CCs including their photo activation are fully determined by a number of external 
parameters, namely (Q,R,e*,T). Before we discuss the dynamics of the CC with the 
initial conditions of the electrons in phase space (pj,fj), we elaborate on a universal 
scaling of driven Coulomb explosion dynamics. It is based on the fact that all potentials 
(including the ionic background potential) originate ultimately from the homogeneous 
Coulomb forces. It will reveal the dynamics is the same for parameter sets emerging 
from an arbitrary reference set (Q,R,e*,T) by a scaling transformation which we will 
specify. 

The system Hamiltonian Eq. (2) exhibits a global scaling with scaled variables x 
according to 

(p,r) = (r^p,^) (E,t) = (rj E , i]~ 3 / 2 1) , (15) 

which applies to Coulomb systems. Note, that the external potential Eq. (1) belongs to 
this class since it is given as the Coulomb potential of an extended charge distribution. 



i The Heaviside step function is 0(f) = I ? * < ^ 

w [ 1 for t > 
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These scaling properties can be used to identify sets of different external parameters 
which will lead to the same driven Coulomb explosion dynamics, (Q,R,e*,T) — > 
{Q,r ] - 1 R,r]£*,r]- 3 / 2 T). E.g., with 77 = (T/T*) 2 / 3 we can map 

(Q,i?,e*,T)-^(Q,r ? - 1 J R,^*,T*), (16) 

where T* is any reasonable chosen reference pulse length. It should be emphasized, that 
the activation Hamiltonian Ha in Eq. (14) has the same scaling properties, since the step 
function 6 does not posess a time scale, i.e., &(t — r) = (r/ 3//2 (i — f)) = &(t — f). 

We illustrate this scaling property in Fig. 2 where we have shown the converged 
electron spectra dP/dE of three CCs with parameters related by the scaling property 
Eq. (16). In Fig. 2a each spectrum displays the same basic features, mainly a narrow peak 
near E w and a broad peak in the negative energy region, resulting from a complex 
dynamics of electron emission and equilibration of the remaining electron plasma in the 
activated CC which we discuss in Sect. 5. The scaling property becomes manifest in 
Fig. 2b, where the spectra are rescaled according to i]dP(i]^ 1 E)/dE to reveal that they 
are indeed identical within numerical accuracy. 




-120-100 -80 -60 -40 -20 20 40 -120-100-80 -60 -40 -20 20 40 

E [eV] E [eV] 



Figure 2. Electron spectra dP/dE for reference parameters (Q, R [ao], £* [eV], T [fs]) = 
(123,26,10,10) (red, full), (123,13,20,3.54) (green, dashed) and (123,52,5,28.28) 
(blue, dotted) . The latter two parameters are related to the former by scaling factors of 
1] = 2 and 77 = 1/2 respectively. Panel (a) shows the actual electron spectra dP/dE and 
panel (b) the rescaled spectra r)dP{r]~ 1 E)/dE that demonstrate the scaling invariance 
of the underlying dynamics to within numerical accuracy. 



5. Time-resolved electron spectrum 

In order to understand the general feature of a double peaked final electron spectrum 
for an activated CC, as shown in the previous section, it is necessary to appreciate that 
CCs are open systems and as such are susceptible to electron loss. Indeed, it is the 
emission of energetic electrons that facilitates the relaxation of the remaining electrons 
from the highly excited activation state. This comes about as the remaining electrons 
in the CC lose the interaction energy with respect to emitted electrons, leading to a 
deeper effective binding potential. 
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Insight into the formation of the final electron spectrum of a CC is obtained by 
tracing its time-resolved evolution from the well specified initial state described above. 
This is shown in Fig. 3 for parameters (Q,R[&o],e* [eV],T[fs]) = (1000,30,50,20). In 
this spectrum all electrons are considered, activated as well as not yet activated in which 
case the designated kinetic energies according to Eq. (12) are already included. Therefore 
the initial spectrum is mono-energetic with all electrons having a single-particle energy 
of e*. 



trap 




■ i . i i . . I u o.l 

-20 20 40 60 80 100 

t[h] 



Figure 3. Time evolution of an activated Coulomb complex with the following 
parameters (Q, R [a ], e* [eV], T [is]) = (10 3 , 30, 50, 10). (a) Electron emission rate 
through direct emission before trapping time itrap (red) and through thermal emission 
from a trapped plasma after i t rap (blue) and activation rate of electrons (green), (b) 
Time-resolved electron spectrum of Coulomb complex including as yet not activated 
electrons, (d) final electron spectrum at t = 100 fs, separated by activation time of 
electrons before ttrap (direct, red) and after itrap (thermal, blue) as well as total 
spectrum (green). Gray shaded area marks the analytical value for a sequential 
spectrum (see Appendix) . (c) Mean energy of the electrons in time from a Maxwell- 
Boltzmann fit to the kinetic energy spectrum of electrons (see (Fig. 4a)) within a sphere 
of radius 2R (purple) and expressed through the energy scale E (see Fig. 4b) from the 
exponential fit to the spectrum of the thermally emitted electrons (blue). 
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Directly emitted electrons 

In the early stage of the activation, with a rate shown in green in Fig. 3a, activated 
electrons have a positive energy and they may leave the CC, a process which we call 
direct electron emission. With each emitted electron the total charge of the CC increases 
thereby deepening the trapping potential for the subsequently activated electrons, which 
is visible as the decrease in energy of the main peak in the time-resolved electron 
spectrum. 

The direct emission process continues until the main peak of mostly still dormant 
electrons falls below the threshold of e = and newly activated electrons find themselves 
bound to the CC, which they can no longer leave. We refer to this as trapping of electrons 
[9-11] by the now positively charged CC. In terms of the direct emission rate, which is 
shown in red in Fig. 3a, this corresponds to a sharp drop-off. We call the time when the 
emission rate reaches zero the trapping time t trap which here occurs near — lOfs. The 
process of direct emission of electrons is identical to the previously described multi-step 
ionization [12]. It results in a plateau-shaped spectrum shown in red in Fig. 3d. The 
sequential process is amenable to an analytic description giving an exact expression 
for the height of the plateau of purely geometrical origin as AR, with A 0.84 (cf. 
Appendix A for details). The number of directly emitted electrons may therefore be 
readily approximated as A^ direct ss Ae*R which here amounts to less than 5% of the total 
number of electrons. 

Plasma formation and equilibration 

With the onset of trapping at t — t tiaiP the activated electrons within the CC undergo 
energy-exchanging collisions that lead to relaxation. This has two consequences: on 
the one hand the spectral peak undergoes a broadening as a thermalized plasma is 
formed. On the other hand particularly fast electrons in the plasma may leave the 
complex contributing to a second emission peak in Fig. 3a and resulting in a continued 
decrease of the mean energy of the plasma. Since emitted electrons carry away energy 
the mean electron energy in the complex is reduced. This opens a gap in the spectrum 
that individual electrons need to overcome in order to leave the cluster. However, the 
plasma temperature eventually does not suffice to produce fast enough electrons and 
the emission decreases. From this time on the spectrum has essentially assumed its final 
shape. The peak at negative energy describes the thermalized plasma and the peak at 
positive energies contains the electrons emitted early during the activation. The gap 
between these two peaks must necessarily span a multiple of the plasma temperature to 
prevent further electron emission. 

The kinetic energy distribution of the plasma electrons follows closely a Maxwell- 
Boltzmann distribution, as can be seen in Fig. 4a. Since the photo activation continues at 
the time for which the distributions are determined, we consider only electrons activated 
before t — tn, where tn = 2ir(R 3 /Q) 1 / 2 is the oscillation period within the Jellium 
potential which provides an intrinsic time-scale of the electron dynamics. Obviously, 
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Figure 4. (a) Kinetic energy spectra of plasma electrons dP(eki n )/cfckin and fitted 
Maxwell-Boltzmann distributions ~ \J E f^ksT 3 exp{~E/k^T) at different times in 
the evolution of the Coulomb complex, (b) Spectrum dPthmmai/dE of thermally 
emitted electrons at various times and exponential fit ~ exp(—E/E). 



this time is sufficient for the equilibration of newly activated electrons. 

The plasma temperature extracted from the Maxwell-Boltzmann fit exhibits a 
sharp rise beginning at time itrap towards a maximum close to the time of maximum 
activation rate (purple, Fig. 3c). The rise reflects the increasing depth of the positive 
background potential allowing to trap faster electrons in a hotter plasma. Eventually, 
the temperature decreases due to evaporative cooling effects and levels off asymptotically 
slightly below 25 eV. Note, that this temperature remains consistently below the 
temperature inferred from the initial kinetic energy of the electrons, Eq. (12), which 
is due to a broadened spatial distribution of the plasma with respect to the ground 
state. 

Emission from the plasma 

As can be seen in Fig. 4b and as has recently been similarly described elsewhere [13, 14], 
the spectrum dP t hermai/dE of the electrons emitted by the thermalized plasma forms an 
exponential distribution characterized by the energy scale E, 

dP thcrmal /dE = AWna^- 1 exp(-E/E) , (17) 

where A^hermai is the number of electrons thermally emitted from the plasma. Although 
originating from the tail of the equilibrated plasma inside the cluster the "temperature" 
E of the thermally emitted electrons (blue, Fig. 3c) is considerably lower than /c#T of 
the trapped plasma (purple, Fig. 3c). One should therefore not take this experimentally 
directly accessible energy E as the temperature of the trapped electron plasma. 

The exponential distribution Eq. (17) has been measured recently [5] in an 
experiment, where ionization of clusters with 90 eV photons from intense pulses supplied 
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by FLASH [1] was studied by electron spectroscopy. The electron spectra of irradiated 
xenon clusters exhibit high-energy tails of the photo lines, which extend further with 
increasing intensity. These blue wings in the spectra can be related to the formation of an 
electron plasma with supra-atomic density, generated due to the large photo-absorption 
cross section of xenon at 90 eV. Theoretically, these electrons could be distinguished 
from those emitted directly at the beginning of the pulse [5] . Note that electron emission 
may also occur for photo activation just below threshold, e* ~0, as measured recently 
at SCSS [2] for neon clusters [15]. Whereas direct emission is not possible, the energy 
of the excited electrons is redistributed through collisions with some electrons being 
emitted. 

6. Summary 

We have discussed a simple and fairly general scenario for the electron dynamics in dense 
matter exposed to short and intense laser pulses of high intensity. It applies to various 
frequencies u of the light up to X-rays. The lower limit for ui is set by the condition 
that single-photon ionization of many atoms dominates which we call photo activation. 
We have provided a simple formalization of the photo activation in time which relates 
the successive photo-ionization events to the time-dependent shape of the laser pulse 
and the total number N of ionized electrons generated. The process leads inevitably 
to a Coulomb complex, a many particle system of ions and electrons whose Coulomb 
interaction dominates all other forces. From this follows a simple scaling property of 
Coulomb systems which relates seemingly disparate experimental szenarios with respect 
to the size of the system (characterized by the radius R of atomic positions), the time 
scale of photo activation (determined by the pulse width T) and the electron energy 
N(uj — E ip ) made available through the activation. The latter depends through N on 
the total energy contained in the laser pulse. 

Furthermore, we have discussed the time-evolution of the electrons during and after 
photo activation including the formation of a plasma and its equilibration in detail. 
Photo activation in Coulomb complexes is not limited to the examples of intense XUV 
and X-ray pulses illuminating clusters of almost solid density as discussed here. The 
very same phenomenon occurs, e.g., in the activation of ultracold plasmas [16] from 
atom clouds held in magneto-optical traps at appropriately adopted time and energy 
scales. 
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Appendix A. Electron spectrum for sequential ionization 

We derive an analytical expression for the electron spectrum in the case of high excess 
energy e* ^> Q/R and long pulses T — > oo. In these limits all electrons leave the CC 
directly upon activation, without the possibility of exchanging energy with other active 
electrons, i. e., the emission occurs sequentially. Furthermore we assume a homogeneous 
electron distribution within r < R of the non-activated electrons, which gives rise to the 
probability distribution in the radial coordinate for the electron to be activated next as 

dP - 3r2 (A 1) 

Released electrons leave a charged CC behind. This remaining charge, denoted by Q, is 
within our approximation also homogeneously distributed. It modifies the asymptotic 
(or measured) energy E of the next to be activated electron. Therefore the energy 
depends on the radial coordinate r, 

with the modification given by the potential in Eq. (1). The dependence on the 
instantaneous charge Q is emphasized by the subscript. This expression allows one 
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Figure Al. (a) Probability distribution in energy dPg/dE from Eq. (A. 3) of next 
activated electron for CC with parameters (Q, R [ao], e* [au]) = (500, 25, 40) and charge 
Q = 0, 50, . . . , 500 on the CC. (b) Analytical electron spectrum dP/dE from Eq. (A.4) 
from integration of probability distributions. 
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to obtain the probability distribution in terms of the energy as 

dPg = dP 

dE dr 




Q Q\\ 2RJ Q 

where the lower line uses the inverse function of Eq. (A. 2). As the radial coordinate is 
restricted to < r < R, so is the energy in Eq. (A.3), namely —(3/2)Q/R < E—e* < 
—Q/R. The reduction of the energy E—e* is bounded by the potential energy at the 
centre and the surface. The probability distribution of Eq. (A.3) is shown for the specific 
case of e* = 40 au, R = 25 ao und various instantaneous charges Q = 0, 50, ... , 500 in 
Fig. Ala. The larger the charge on the CC is, the wider the distribution in energy 
becomes and the further it moves to lower energy. 

As can be seen in Fig. Al only certain values of Q contribute to the spectrum at 
a given energy E. These values are given by Q\{E) < Q < Q2(E) with Q\{E) = 
(2/3)R(e*—E) and Q2{E) = min (R(e*— E), Q). Remember that Q is the total charge 
reached after the activation and removal of all electrons. Thus one obtains by integration 

*P = [® 2{E) d Q dP Q 
dE Jqi(e) dE 

A R for -Q/R < E-e* < 0, 

\q,r{E) for -(3/2)Q/R < E-e* < -Q/R. { ' } 

Interestingly the electron spectrum for higher energies is independent of energy E 
with a system-independent constant given by A = 3 [2^3 In (3+v^) -y/3hx (6) -2] w 
0.843. Towards lower energies E it falls off monotonically according to Xq^(E) = 
3R [2V3\n (3+v^6x) -y/3 In (6x) -2v^2^] with X = (e*-E)R/Q. 

The entire analytical spectrum is shown in Fig. Alb for the specific parameters 
Q = 500, R = 25 ao und e* = 40 au. In particular the plateau in the region from 




[eV] 



Figure A2. Comparison of numerical (symbols) and analytical (lines and shaded) 
electron spectra for various parameters. Analytical spectra calculated according to 
Eq. (A. 4). A pulse length of T = 10 fs, sufficiently long for sequential ionization, was 
used for all three spectra. 
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s* — Q/R to e* is a characteristic feature of the sequential ionization as seen in an 
experiment [12] and numerical simulations [17]. While the height of the spectrum is 
solely determined by the radius R, the width of the full spectrum and of the plateau- 
region are given by 3Q/2R and Q/R respectively. Note that this implies that in the 
case of e* < 3Q/2R trapping of electrons will occur and the full sequential spectrum 
is not realized. However, for those electrons emitted before the onset of trapping the 
sequential spectrum as derived above is valid. 

We compare the analytical spectrum of Eq. (A. 4) with numerical results for three 
sets of parameters in Fig. A2, in each case with a pulse length of T = 10 fs which proves 
sufficiently long for the CC to exhibit sequential ionization behaviour. Indeed, we find an 
overall excellent agreement between the numerical and the analytical result, in particular 
as pertains to the appearance of an extended plateau region, well described in width, 
height and position by Eq. (A. 4). Nevertheless, some minor discrepancies owing to the 
continuous approach employed in the analytical calculation can be observed. A narrow 
peak at e* in the numerical spectrum arises from the first activated electron with exactly 
the excess energy e*. Furthermore, in the energy domain —3Q/2R < E—e* < —Q/R, 
the shell-like structure of the initial configuration of the CC, as shown in Fig. 1, leads 
to a step-like spectrum, as opposed to the smooth increase in the analytical case. 
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